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ABSTRACT 

The early validation of a 20 year service integrity for the bonded 
interface in solar cell modules is an important requirement in the Low Cost 
Solar Array (LSA) project. 

The first annual report (Science Center Report No, SC5106.22AR) out- 
lines and implements a physical /chemical evaluation program for solar cell 
encapsul ants. The results of computer controlled ultrasonic and optical /el 1 ipso- 
metric mapping for interfacial defect characterization in solar modules is 
summarized in the second annual report {SC5106.49AR). The developnent and 
validation of an atmospheric corrosion model and test pin for LSA outdoor 
service at the Mead, Nebraska test site is presented in the third annual 
report (SC5106. 86AR) . 

In the present fourth phase of study detailed in this annual report 
emphasis is placed on the development of AC impedance as a nondestructive 
evaluation (NDEl methodology for solar arrays and the further development of 
corrosion models and materials selection criteria for corrosion resistant 
interfaces. 
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1.0 OBJECTIVES 


The general objective of this phase of the program is to broaden the 
present corrosion model and atmospheric corrosion studios developed for the 
Mead, Nebraska test site. The two major objectives Included in this present 
study are stated as follows: 

1) Continue further development of atmospheric corrosion monitors and AC 
impedance measurements as nondestructive evaluation (NDE) tools for 
LSA module perfonriance and life prediction. 

2) Initiate development of materials selection criteria and tests for 
corrosion protection and environment resistant interfaces as required 
for validation of corrosion models. 
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2»0 SUMMARY 

Atmospheric Corrosion Model Verification 

a. A 13 month experimental study of climatology and atmospheric 
corrosion monitoring (ACM) of the Mead, Nebraska ISA test site 
has been compl .-ted. 

b. During this 13 month period the ACM units recorded the corrosion 
protection function of an encapsulant system (a reactive primer 
GE-SS4155 on Zn/Cu plates encapul ated by Sylgard 184). 

c. A non- encapsul ated ACM unit correlated Mead climatology with 
corrosion rates as defined by an atmospheric corrosion model. 

d. Both the intermediate and final results of this study verify the 
prediction that atmospheric corrosion rate is the product of the 
moisture condensation probability (P^,) and the maximun ionic 
diffusion current (l^) at the corrosion surface or interface. 

e. Encapsulant corrosion protection is specifically related to the 
suppression of l|^ at the potential corrosion interface. 
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2. 2 AC Impedance Monitoring and Analysis 

a. A predictive model for solar array current- voltage (I-V) 
response has been developed and correlated to AC impedance. 

b. A new and simplified methodology for conducting AC impedance 
measurement and spectruii analysis has been developed and 
verified for nondestructive evaluation of LSA performance. 

c. A preliminary design for an AC impedance measurement, spectrua 
analysis, and performance optimizing control for a solar array 
branch circuit has been outlined. 

H. 3 Hydrothermal Stress Analysis of Solar Arrays 

a. A computer model for hydrothermal stress analysis (HTSA) has 
been successfully applied to evaluate the combined effects of 
temperature- humidity cycling on solar cell cracking mechanism in 
LSA modules using a fiber board substrate. 

b. The protective effect of moisture diffusion barrier coatings of 
ethylene vinylacetate (EVA) and polytrifluorochloroethylene 
(KEL-F) is predicted to delay but not change the solar cell 
cracking process. 
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Criteria for Encapsulant Bonding 

a. A new criterion for encapsulant bonding is being developed and 
specifically directed at supressing the micrccorrosion process 
in the presence of internal defects such as micro- cracks. 
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3.0 INTRODUCTION 


In view of the promising early results of the atmospheric corrosion 
model developed for the low cost solar array (LSA) at the Mead, Nebraska test 
site the present program was expanded relative to a general LSA life predic- 
tion program as shown in Fig. 1. Tlie logic flow diagram of Fig. I describes 
the materials imputs to an LSA system in the upper block. The next lower 
block describes the environmental stresses known to cause LSA performance 
degradation. The specific mechanisms of degradation and consequent degrada- 
tion effects ere detailed in the subsequent decending blocks of Fig. 1 with 
debonding and corrosion near the bottan of the figure and directly connected 
with photovoltaic reliability and durability. On the right of Fig. 1 are 
categorized the specialized measurementand analysis programs which correlate 
v/ith each degradation mechanism and effect. 


The prior armual report (505106. 86AR) details the results of the 
corrosion program whose scope is defined by the inner bracket on the right of 
Fig. 1. In this prior program consistant interrelations have identified the 
corrosion chemistry, changes in the AC impedance spectrum, and the photo- 
voltaic current- vol tage (IV) profile which measures photovolaic performance. 

In the present program the corrosion studies were expanded to include develop- 
ment of a life prediction model and life prediction test plan which included 
the debonding and corrosion process. Hi is report details the progress of this 
effort under a Task I which now encompasses corrosion modeling, life predic- 
tion in terms of photovoltaic performance and nondestructive evaluation of 
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Fig. 1 General LSA life prediction program. 
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performance of modules and solar array branch circuits. This phase of the 
study is closely coupled to extensive studies in photovoltaic design optimi- 
zation already completed by JPL. 

In order to provide for new materials selections in terms of 
advancing LSA technology a new study program in material selection criteria 
for corrosion resistant bonding in the presences of microdefects has been 
initiated. This study is defined as Task 2 in the right bracket of Fig. 1. 
This task Implements an already developed computer mcdel^^^ for hydrothermal 
stress analysis (HTSA) to define potential mechanisms of LSA microfailure and 
to focus the materials studies and bonding criteria to the specific failure 
mechanisms. In addition, a new physical/chemical borviing criteria and test 
protocol has been initiated for minimizing ionic conduction at encapsulant 
interfaces in LSA. 
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4.0 ANALYTICAL MODELS 

Twd analytical models have been developed v/ithin this program to 
translate the results of corrosion analysis into predictions of photovoltaic 
performance and durability. The first of these models develops the current- 
voltage relations for a single solar cell performance. By introduction of 
summing relations for series and parallel interconnection the performance and 
AC impedance response of a branch circuit can be calculated. Tlie second model 
is developed for a simplified computational methodology to analyse AC imped- 
ance spectron. The combination of these models provides a potential new basis 
for nondestructive evaluation (NDE) and feedback control of solar arrays. 

4. 1 Sol ar Array Current Voltage (I-V) Model 

Appendix I provides a detailed derivation of an analog model for 
solar cell response which permits direct calculation of the I-V fill factor, 
fiiaxlinum power point and other performance properties relating to solar cell 
response. The equation, diagram and nomenclature of Fig. 2 give a brief 
synopsis of this I-V response model. Based upon controlled variations of 
Isolation S and temperature T the curves of Fig. 3 show how the design para- 
meters of this I-V response model can be evaluated experimentally. The I-V 
curves in the upper right portion of Fig. 3 are taken from a published report 
and represent typical silicon module response. The calibration experiments 
of Fig. 3 define all of the solar cell design parameters with the exception of 
and Rg. As documented in the third annual report (SC5106.86AR) AC-impedance 
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is very sensitive to small decreases In shunt resistance Rg|^ due to cell 
corrosion which are not evident in the shape of the I-V performance curve as 
defined by the fill factor F. The calculated I-V response curves of Fig. 4 
illustrate why this experimental result is observed. The right view of Fig. 4 
plots the logarithm of versus the calculated fill factor F. For high 

performance cells with F>0.8 the right curve of Fig. 4 shows that F is quite 
insensitive to small changes in R5}^/R5>100. The upper left I-V curve of 
Fig. 4 for KVn=lCA/0 represents a high performance cell response with F“0. 9. 

The slope of the upper curve at V/Vq= 0, which is a measure of -1 /(Rsh+Rs^* 
so slight as to make shunt resistance measurements impractical from I-V 
curves. 

The mathematical selections and schematic curves of Fig. 5 introduce 
an AC impedance model for a single solar cell. Extensive experiments now show 
that single solar cells display a single relaxation time t and when the 
storage V and loss Z" components of the complex impedance are plotted as in 
the lower right curve of Fig. 4a semicircle with a center at Z"=0, Z'= 

Rs + radius Rsh/ 2 results which is tenned the Cole-Cole plot.^®^ 

Evaluation of AC impedance daca via the Cole-Cole plot isolates the shunt and 
series resistance for the single solar cell or a solar cell module as prior 
reports in this program have documented. 

In general, as analytical modeling and direct testing have indicated, 
the presence of mismatches in single cell 1-V response within a solar array 
leads to Towered power efficiency in a solar array. The mathematical 

relations of Fig. 6 illustrate a method for describing statistical varia- 
bilities of all solar cell design properties in terms of Weibull 
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Fig. 4 Calculated current - voltage (I-V) response. 
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cimulative probability functions. The property mean value and distribution 
shape fsctor m^ can be evaluated from ineasuranents or theoretically assigned. 
This process for describing the statistical variability for a single cell can 
be repeated for each of the N cells making up the model branch circuit. The 
I-V response and AC Impedance of each model cell can be calculated and 
combined In appropriate summation relations which describe the series- parallel 
cell interconnections In the branch circuit. Graphical Illustrations for I-V 
summing in series-parallel networks are Illustrated in Fig. 7. The methods by 
which particular types of defects can be located by I-V nieasurement of 
separated portions of a branch circuit are graphically illustrated in Fig. 8. 


Inpedance Model 


The second portion of the analytical modeling effort has been 
directed at developing AC impedance as a nondestructive evaluation (NDE) tool 
for branch circuit analysis as a substitute for I-V characterization. 

Appendix II provides a detailed theoretical discussion and experimental 
appraisal of electrical impedance measurements and impedance spectum analysis 
of solar cells and elements of branch circuits. As pointed out in Appendix 
II, the new development of a simplified least squares circle fit of AC imped- 
ance data to a Cole-Cole distribution function provides an important 
augmentation of the methodology. This circular least-squares method is 
graphically illustrated in Fig. 9. A minimun of three AC data points with 


16 

C3047A/ES 



WEIBULL DISTRIBUTIONS 

p. = CUMULATIVE PROBABILITY INVERTED FUNCTIONS SC81-11663 



17 



-ts* 




Rockwell International 

Sci«nc6 Canter 


SC5I06.104AR 






Fig. 7 I-V curve calculation for series parallel networks (from Ref. 6). 



FAILED CELL' 
(OPEN CELL) 





Rockwell International 

Science Center 


SC5106.104AR 





00 


STORAGE IMPEDANCE Z' 


Fig. 9 Schematic representation of least squares circle fit of the Cole-Cole 
AC impedance fund ton. 
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sufficient seperation on the Z" vs Z* cooixlinates of Fig, 9 are sufficient to 
detennine the values of R 5 and and the distribution quality factor 0 <B< 1,0 
for multiple interconnected solar cells. If all the cells have identical 
impedance response then 3*1,0. As Appendix II shows, the mismatch of AC 
impedance response tends to lower 3 and at the same time produce mismatches in 
I-V responses between cells which lower photovoltaic performance. 


The potential simplicity of AC impedance measurement and spectruii 
analysis as outlined in Appendix II suggest the potential of re<note impedance, 
monitoring and feedback control of a solar array branch circuit. A prel imi- 
nary model for AC impedance measurement, impedance spectrun analysis, and 
feedback control to an LSA branch circuit is graphically outlined in 
Fig, 10. In brief, the- solar array branch circuit is viewed as equipped wiih 
digitally controlled series-parall el interconnects betvveen modules. AC 
impedance measurenients conducted at the solar array site would be digitized 
and tranirnitted to a control location of impedance spectrum analysis. An 
analytical branch circuit IV (current-vol tage) model and predicted impedance 
spectruii wuld be incorporated in a corrective action analysis. Corrective 
action could be transmitted to the power control (PC) unit at the solar array 
site by digital encoded instructions that modify the series to parallel (S/P) 
interconnects in the solar array. The modified S/P interconnection can be 
utilized to locate deteriorated modules and to reoptimize the branch circuit 
electrical outputs for maximun power output. 
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AC impedance Is viewed as a method for augmenting power output 
measurement to provide highly sensitive and interpretable data on environ- 
mental aging effects. Tbe interpretation and control methodology outlined 
above would integrate nondestructive evaluation (NOE) with a corrective action 
program to locate and minimize environmental aging effects on array 
performance. 
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5.0 EXPERIMENTAL 

This section discusses and interrelates the results of natural corro- 
sion studies at the Mead test site and computer modeling and accelerated aging 
at the Science Center. 

5. 1 Mead Site Atmospheric Corrosion Monitoring 

As reported in the third annual report (SC5106.86AR), a new atmos- 
pheric corrosion model, as outlined in Fig. 11 was developed and implemented 
in corrosion studies at the Mead test site. As shown in Fig. 11, the corro- 
sion rate is predicted to be the product of the surface condensation proba- 
bility of water vapor Pq and the diffusion controlled corrosion current 
Thi’s model combines surface physical chefnistry and electrochemisty and is 
verified by the direct correlation between relative humidity and the logarithm 
of measured corrosion current as shown in the upper right view of Fig. 11. 

The corrosion monitor experiment was initiated at Mead, Nebraska on 
July 12, 1979 and continued for over a year until terminated in October 1980. 

A pnotographic view of the Mead, Nebraska LSA test site is shown in the upper 
view of Fig. 12. The lower view of Fig. 12 shows two atmospheric corrosion 
monitors (ACM) installed on Ine right rear portion of the array. The monitors 
consist of the edges of alternating plates of Zn and Cu seperated by a 50 m 
polyester (DuPont-Myl ar) dielectric film. Moisture condensation on the 
exposed surface closes the circuit between the plates and current flows. 

The solid state circuitry of the ACM converts the logarithm of the current 
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Fig. 12 Photographic view of Mead, Nebraska LSA test 
site (upper view) and atmospheric corrosion 
monitors (low^r view) installed on rear right 
portion of array. 
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(microamperes) to a proportional DC voltage which is transmicted to the 
recording station located in the power control unit (PCU) located behind the 
right end of the array (see upper view Fig. 12). Specific details of the 
corrosion and climatology recordings are discussed in the third annual report 
(SC5106.86AR). 

During the 13 month deployment at the Mead, Nebraska test site the 
ACM units recorded, at 10 minute intervals, the corrosion protection function 
of an encapsulant system consisting of a reactive primer. General Electric 
GE-SS4155, on Zn/Cu plates encapsulated by 2nin of Dow Sylgard 184 silicone. A 
non-encapsul ated ACM unit correlated Mead climatology with corrosion rates as 
defined by the atmospheric corrosion model outlined in Fig. 11. 

The fundamental assumptions of the corrosion model appear verified in 
this 13 month study. As shown in the upper curves Fig. 13 and Fig. 14 the 
encapsulated ACM showed a negligible corrosion rate under both wet and dry 
conditions during the entire 13 month exposure period. The nonencapsul ated 
ACM unit showed a linear interrelation between the logarithm of corrosion 
current I and moisture supersaturation temperature (Tq-T) (dew point minus 
temperature) as defined by the corrosion model. Below a critical relative 
humidity or (Tq-T) for surface condensation the condensation probability is 
negligible and corrosion current is minimized to the lower limits of current 
detection. During periods of surface condensation the uncoated ACM achieves a 
maximum current output which is the diffusion controlled current I|_. The 
negligible current of the encapsulated ACM identifies the mechanism of 
encapsulant corrosion protection as specifically related to the suppression of 
II at the potential corroding interface. 
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Outdoor aging effects on the inaxiinun corrosion current at lOOX R. H. 
of the uncoated ACM are shown in Fig. 15. This maxiniun current displays 
weekly and seasonal variability which depends upon the ionic conductivity of 
surface film produced by both corrosion and dirt accumulation. The upper 
curve of Fig. 15 points out that this maximum current decreases from « 20 uA 
for an unaged ACM to «5 uA (".fter 13 months exposure. As previously 
docunented and shown in Fig. 15, the freeze- thaw cycles of Mead winter 
cause a wide variability in maximiin current due to lowered diffusion and ion 
solubility in ice as opposed to water. 

The curves of Fig. 16 clearly show the effects of Mead aging by 
plotting the logarithm of corrosion current I versus moisture supersaturation 
temperature (Tq - T). These plots are obtained by parametric cross plotting 
the respective I and (Tq - T) data of Fig. 13 and Fig. 14 at equivalent time 
conditions. The upper curve of Fig. 16 shows a lower limit corrosion current 
I «0.04 uA characteristic of the lower limit of measurement sensitivity. At 
this lower limiting current the dry uncoated and encapsulated ACM are equiva- 
lent in response. As shown in Fig. 16 with increasing supersaturation temper- 
atures (Tq - T) the corrosion current displaces to an upper limit character- 
istic of diffusion controlled response. Referring to the corrosion model (see 
Fig. 11) the transition region for corrosion current output is related to 
water wettability and capillary absorption. The curves of Fig. 16 show that 
this transition region is significantly shifted to higher values of Tq - T as 
a consequency of t4ead outdoor exposure. The same aging factors which tend to 
lower maximun current, such as film formation and contamination, can be 
expected to diminish water wettability and capillarity. 
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CUMULATIVE PAYS MEAD OUTDOOR EXPOSURE 



DAY OF YEAR 

YEAR 1979 YEAR 1980 


Fig. 15 


Maximum atmospheric corrosion monitor (ACM) current I at 100% 
relative humidity versus Mead site exposure time (days). 
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(SEPT. 14-16, 1980) 
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Fig. 16 Corrosion monitor current (I) vs. moisture super saturation 

temperature (Tp-T) at 28 days (upper view) and 429 days (lower view) 
outdoor Mead exposure. 
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Photographic views of the uncoated and encapsulated ACM units 
subsequent to 430 days Mead outdoor exposure are shown in Fig. 17. The 
uncoated ACM shows ample evidence of extensive corrosion and material erosion 
while the encapsulated ACM shows the original metallic luster and intact 
interfacial bond to encapsul ant characteristic of complete corrosion suppres- 
sion at the Cu/Zn plate surface. Based upon these results the 6E-SS4155 
reactive primer appears as a strong condidate in planned studies of interface 
corrosion durability. 

5. 2 Hydrothermal Stress Analysis 


n 2l 

An earlier environmental study by E. Laue and A. Gupta at JPL' ‘ 
showed that combined moisture and UV cycling of an EVA (ethylene vinylacetate) 
encapsulated silicon solar cell on a wood composite substrate produced cell 
cracking and evidence of subsequent microcrack corrosion. The Science Center 
had recently completed development of a hyrothermal stress analysis (HTSA) 
program for composite laminates^^^ and a study was undertaken to explore the 
applications of the HTSA model in ISA life prediction. The logic flow chart 
for the HTSA model is shown in Fig. 18. As shown in Fig. 18 the HTSA model 
imputs a physical description of layers, the hydrothermal , and mechanical 
properties of each layer, and the discription of the hydrothermal cycles 
imposed in the simulated aging study. 

A cross section of the solar cell laminate considered for HTSA 
modeling is shown in the upper view of Fig. 19. Where the cell front surface 
faces right. The combined variations in temperature and humidity are shown in 
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Fig. 17 Uncoated (left) and er»capsu1ated (right) corrosion monitors 
after 430 days outdoor exposure at Mead, Nebraska test site. 
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A. DENSITY, EXP. COEF.aj.an 

B. THERMAL CONDUCTIVITY 

C. DIFFUSION CCEF. (TEMP. DEP.) 

D. INITIAL CONDITIONS 
SURFACE TEMP. & R.H FOR EACH 
TIME CYCLES 

number OF TIME CYCLES 



PRINTOUT OF INPUT 
VARIABLES 
FORAC110 / 




INPUT AC3 VARIABLES 
CALCULATE LOCATION 
OF MIDDLE LAYER 
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18 Flowchart of hydrothermal stress analysis program 












Rockwell International 

Sclance Center 


SC5106.104AR 


INTERVAL 



HOURS 


Fig. 19 Hydrothermal history of one environmental exposure cycle. 
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the lower curves of Fig. 19, The first computer run did not include layer 5 
which represents a moisture diffusion barrier film on the back surface of the 
solar cell. 

Moisture diffusion properties for both clear and pigmented EVA 
encapsul ant films obtained from JPL were experimentally detennined and these 
properties are reported in Table 1, The full temperature dependent moisture 
and thennal diffusion imputs for each of the four solar cell layers (see 
Fig. 19) which are input to the HTSA model are summarized in Table 2, The 
hydro- thennal- elastic coefficients for the four solar cell layers for HTSA 
input are detailed in Table 3. 

The computed moisture and stress profiles produced by the end of 
first 48 hours interval of hydration for the four layer laminate is shown 
upper Fig. 20. The moisture is seen to penetrate both layer 1 and 4 and to 
produce a substantial biaxial tensile stress in the silicon wafer which is 
counterbalanced by a lower biaxial compressive stress in the fiber board. 
Encapsul ant layers 1 and 3 show minor biaxial tension at this point. 

The lower view of Fig. 20 shows the moisture and stress profiles at 
the end of the 24 hr drying cycle which fonns the second interval of aging. 

All the moisture in layer 1 is essentially removed but only the surface mois- 
ture is lowered in the thicker fiber board layer 4. Both the tensile stresses 
in the silicon layer 2 and compressive stress in layer 4 are slightly dimi- 
nished by this drying cycle. An interesting and important feature of inter- 
laminal cracking in the fiber board layer 4 is shown 1n the steep moisture 
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Table 1 Moisture Diffusion Properties of Solar Cell Encapsulants 
Diffusion Equation 0 ■ Dq exp (-Eq/RT) 


Moisture diffusion coeff. 
at 22“C, (mr/sec) 

5 X 10“^ 

1 X lO'^ 

Max. HoO uptake at 100% RH 
R~T 

0. 08% 

0,3% 

Activation energy, Eq, 
kcal/mole 

12 

13 

Dq (nil? /sec) 

316 

370 
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Table 2 Moisture Diffusion Input to HTSA Model 

‘!*JHBERS OF LAYERS * 4 

THE PLATE IS EXPOSED ON BOTH SIDES. 


INITIAL COCENTRATION DISTRIBUTION 
LAYER I C(0/0) AT X(IN) 

X .0000 .0045 .0095 .0143 ,0190 

C .0000 ,0000 .0000 ,0000 .0000 

LAYER 2 C(0/0) AT X(1N) 

X ,0000 .0047 .0095 .0143 .0190 

C .0000 .0000 .0000 .0000 .0000 

LAYER 3 C(0/0) AT X{IN) 

X .0000 .0047 .0095 .0143 .0190 

C .0000 .0000 ,0000 ,0000 .0000 

LAYER 4 C(0/0) AT X(IN) 

nnnn 

• iiUUU .U4U ♦OVJO eliiOO 

.0000 .0000 .0000 .0000 .0000 

LAYER 1 

LAYER THICKNESS H • .0190 INCHES 
THERMAL CONDUCTIVITY • .15000 BTU/(m-FT-F) 

DENSITY . 57.00000L8M/(FT**3) 

FOR TEMPERATURE LESS THAN .0 DEGREE C 
FOR TEMPERATURE GREATER THAN .0 DEGREE C 
TEMP. IS IN DEGREE K 

BOUNDARY CONCENTRATION IS GIVEN BY EQUATION 
LAYER 2 

LAYER THICKNESS H - .0190 INCHES 
THERMAL CONDUCTIVITY » 1.00000 BTU/OR-FT-F) 
DENSITY » 144.00000LBM/(FT**3) 

FOR TEMPERATURE LESS THAN .0 DEGREE C 
FOR TEMPERATURE GREATER THAN .0 DEGREE C 
TE!1P. IS IN DEGREE K 

BOUNDARY CONCENTRATION IS GIVEN BY EQUATION 
LAYER 3 

LAYER THICKNESS H » ,0190 INCHES 
THERMAL CONDUCTIVITY - .15000 BTU/OR-FT-F) 

DENSITY » 57.00000LBM/FT**3) 

FOR TEMPERATURE LESS THAN .0 DEGREE C 

FOR TEMPERATURE GREATER THAN .0 DEGREE C 
TEMP. IS IN DEGREE K 

BOUNDARY CONCENTRATION IS GIVEN BY EQUATION 
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(4 Layers, Fig. 19) 


D - ,2Q0E-02EXP( -3000. 0/TEHP. ) MM^*2/SEC. 

D - .570E OOEXP( -5000.0/TEMP.) MM**2/SEC 

C » .001500 mUM.** 1.00 

0 » .100E-04EXP( -500000.0/TEMP.) MM**2/SEC. 

0 » .570E OOEXPC -5000,0/TEMP.) MH**2/SEC, 

C » .00000 NIUM.** 1.00 

D » .200E-02EXP( -3000. 0/TEMP. ) W4**2/SEC. 

D • .570E 00EXP( -5000,,0/TEMP.) Mm*2/SEC. 

C - .003000 mUH.** 1.00 
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Table 2 (continued) 


LAYER 4 

LAYER THICKNESS H • .1250 INCHES 
THERMAL CONDUCTIVITY » .33000 3TU/(IR-FT-F) 

DENSITY » 72.00000LBM/(FT**3) 

FOR TEMPERATURE LESS THAN ,0 DEGREE C 0 * .200E-01EXP( -L6000. 0/TEMP.) MM**2/SEC. 

FOR TEMPERATURE GREATER THAN .0 DEGREE C D • .570E 00EXP( -500D.0/TEMP.) MM**2/SEC. 

TEMP. IS IN DEGREE K 


BOUNDARY CONCENTRATION IS GIVEN BY EQUATION C » .250000 *HLW.** l.OO 


TIME 

START 

START 

FINAL 

FINAL 

START 

START 

FINAL 

FINAL 

HOUR 

LEFT T 

RIGHT T 

LEFT r 

RIGHT T 

LEFT H 

RIGHT H 

LEFT II 

RIGHT H 

.0 TO 48,0 

212.0 

212.0 

80,0 

80.0 

.0 

.0 

100.0 

100.0 

48.0 TO 72.0 

80.0 

80.0 

212. 0 

212.0 

100.0 

100.0 

.0 

.0 

72.0 TO 120.0 

212.0 

212.0 

80.0 

80.0 

.0 

.0 

100.0 

100.0 

120,0 TO 144.0 

80.0 

80.0 

212.0 

212.0 

100. 0 

100.0 

XXX. 0 

XXX. 0 


NUMBER OF CYCLES 10 
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Table 3 Hydro-Thermal-Elastlc Imput to HTSA Model (4 Layers, Fig. 19) 
POINT STRESS ANALYSIS 

(REFERENCE CH. 2.1 OF THE ADVANCED COMPOSITE DESIGN GUIDE - THIRD ED.) 


INPUT DATA 

NL • 4 MM » 4 NS - 1 

PRINT DATA FOR UYERS I, 4, I 
THE STRESS FAILURE CRITERIA IS USED 

MATERIAL 1 IS CLEAR EVA 

E(L) » 20000. E(T) - 20000. G(LT) -» 

U(LT) ■ .5000 ALPHA(l) » .0002000 ALPHA{T) - 
DEFAULTS TAU «< .0 THICKNESS - 


NOMENCLATURE 

NLD - 1 NL « NO. LAYERS 

M4 - NO. MATERIALS 
NS » NO. SETS 
NLD - NO. LOAD CONDITIONS 
E « YOUNGS MODULUS 
U - POISSONS RATIO 

600. ALPHW(L) » .0000800 ALPHA-COEF THERMAL EXPANSION 

.0002000 ALPHW(T) » .0000800 EPS=AXIAL STRAIN 

.0200 CHI -ROTATIONAL STRAIN 


SIGA(L) « 1890. SIGA(T) " 1890. SIGA(LT) » 110. 

SIGA(L) - -1890. SIGA(T) » -1890. N,M,Q-AXIAL, MOMENT, SHEAR LOAD 


MATERIAL 2 IS SILICON 


E(L) « 24000000. 

E T) « 

24000000. 

G(LT) » 10000000. 

ALPh-W(L) 

sc 

.0000010 


U(LT) » .2100 

ALPHA(L) * 

.0000020 ALPHA(T) - .0000020 

ALPHW(T) 

* 

.0000010 


DEFAULTS 

TAU - 

.0 THICKNESS » .0200 




SIGA(L) » 40000. 

SIGA(T) 

» 40000. 

SIGA(LT) « 10000. 





SIGA(L) - -40000. 

SIGA(T) 

» -40000. 




N,M,Q»AXIAL, MOMENT, SHEAR 

LOAD 

MATERIAL 3 IS 

EVA WHITE 






E(L) * 20000. 

E(T) = 

20000. 

G(LT) * 600 
ALPHA{T) = .0002000 

ALPHW(L) 

a 

. 0000800 


U(LT) - .5000 ALPHA )L » 

.0002000 

ALPHW(T) 

s 

. 0000800 


DEFAULTS 

TAU - 

.0 

THICKNESS » .0200 




SIGA(L) - 1710. 

SIGA(T) 

» 1710. 

SIGA(LT) - no. 





SIGA(L) » -1710. 

SIGA(T) 

» -1710. 



N,M,Q=AXIAL, MOMENT, SHEAR 

LOAD 

MATERIAL 4 IS FIBERBOARO 






E(L) 5 900000. 

E(T) - 

900000. 

G(LT) = 300000. 

ALPHW(L) 

a 

.0063000 


U(LT) * .3000 ALPHA(L) - 

.0000040 

G(LT) » .0000040 

ALPHW(L) 

8 

.0063000 


DEFAULTS 

TAU -> 

.0 

THICKNESS =• .1250 




SIGA(L) = 4430. 

SIGAIT 

) a430. 

SIGA(LT) = 2860. 



N,M,Q-AXIAL, MOMENT, SHEAR 

LOAD 

SIGA(L) =■ -4620. 

SIGA(T 

) - -4620. 





NT,MT=*N,M*DUE TO TAU 

TAU-TEMPERATURE DIFFERENTIAL 
SIG =• APPLIED STRESS 
EPSA = ALLOWABLE STRAIN 
SIGA = ALLOWABLE STRESS 
SIGBAR = AVERAGE STRESS 
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concentration gradient near the surface in lower Fig. 20. The important cunu- 
lative effect of repeating these hydro-thermal exposure cycles is shown in the 
curves of Fig. 21. 

Repetition of the hydrothermal cycles shown in Fig. 19 leads to a 
cumulative buildup of moisture concentration in the fiber board layer 4. The 
resultant increased biaxial compressive stress in layer 4 is counter balanced 
by Increasing biaxial tension in the silicon solar cell in layer 2. The life 
prediction presented by the HTSA model is evaluated by the stress a in each 
layer divided by the allowed stress o( All owed) which represents the strength 
reduced by a margin of stress representing the allowed failure probability. 

The upper curves of Fig. 21 show the safety margin for the silicon wafer 
diminish as moisture accumulates in the fiber board in the lower curves of 
this illustration. The four layer solar cell is predicted to fail in about 
three six day cycles of alternate hydration and drying. 

The effect of a fifth layer cell layer of a thin (0.003 in) moisture 
barrier coating on the predicted moisture build and safety margin curves are 
shown in the two right curves of Fig. 21. The HTSA inputs which describe 
these protective coatings are described in Table 4. Addition of this thin 
layer of EVA or the less permeable Kel-F coating is shown in Fig. 21 to sub- 
stantially delay but not change the failure mechanism described for the 
uncoated solar cell. The Kel-F barrier coating was chosen for its extremely 
low moisture permeability coefficient which is 300 times lower than polyethylene. 
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CYCLES (1 CYCLE -6 DAYS) 



F1g. 21 Moisture content and safety margins of fiberboand in encapsulated 
sol ar cel 1 . 
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Table 4 HTSA Model Iinputs for Moisture Barrier Coatings (Layer 5, Fig. 19) 

LAYER 5 - CLEAR EVA 

LAYER THICKNESS H • .0030 INCHES 

THERMAL CONDUCTIVITY « .15000 BTU/(m-FT-F) 

DENSITY - 57.00000L0M/(FT**3) 

FOR TEMPERATURE LESS THAN .0 DEGREE C 0 - .200E-02EXP( -3000.0/TEMP.) MM**2/SEC. 

FOR TEMPERATURE GREATER THAN .0 DEGREE C D - .570E 00EXP( -5000.0/TEMP.) MM**2/SEC. 

TEMP. IS IN DEGREE K 

BOUNDARY CONCENTRATION IS GIVEN BY EQUATION C - .001500 *HUM.** l.OO 
MATERIAL 5 IS CLEAR EVA 


E(L) • 

20000. 

E T) - 
ALPHA (L) - 

20000. 

G(LT) - 600. 

ALPHW(L) • .0000800 

U(LT) - 

.5000 

.0002000 

ALPHA(T) ■ .0002000 

ALPHW(T) ■ .0000800 

DEFAULTS 

TAU - 

0 

THICKNESS » .0030 

SIGA(L) 

» 1890. 

SIGA(T) 

• 1890. 

SIGA(LT) - UO. 


SIGA(L) 

» -1890. 

SIGA(T) 

^ -1890. 

N,H.Q-AXIAL, MOMENT, SHEAR LOAD 


LAYER 5 * KEL-F 

LAYER THICKNESS H « .0030 INCHES 
THERMAL CONDUCTIVITY • .15000 BTU/{m-FT-F) 

DENSITY - 57.00000LBM/(FT^3) 

FOR TEMPERATURE LESS THAN .0 DEGREE C 0- .600E-06EXP( -30000.0/TEMP.) m**2/SEC. 

FOR TEMPERATURE GREATER THAN .0 DEGREE C 0 » .S70E 00EXP( -5000.0/TEMP.) MM**2/SEC. 

TEMP. IS IN DEGREE K 

BOUNDARY CONCENTRATION IS GIVEN BY EQUATION C • .000400 *HUM.** 1.00 
MATERIAL 5 IS POLYTRIFLUOROCHLOROETHYL 


E(L) « 20000. 

E(T) - 20000. 

G(LT) - 

600 

ALPHWIL 

) - .0000800 

U(LT) - .5000 

ALPHA(L) • .0002000 

ALPHA(T) » 

. 0002000 

ALPHW(T 

) « .0000800 

OEFALULTS 

TAU • 0 

THICKNESS - 

.0030 



SIGA(L) - 1890. 

SIGA(T) « 1890. 

SIGA(LT) - 

110. 



SIGA{L) - -1890. 

SIGA(T) - -1890. 



N,M,Q-AXIAL, MOMENT. SHEAR LOAD 
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5. 3 Hydro- thermal Cell Cracking and Corrosion 

Experimental solar cells of the type studied by E. Laue and A. Gupta 
at were provided for study at the Science Center. One cell was hydro- 

thermal ly aged at combined high humidity (1005^ R.H. ) and temperature (75%) 
and cell cracking as predicted by the HTSA model calculations occurred in 8 
days. The type of cells and the pattern of very open cracks in the cell 
surface are shown in the photo views of Fig. 22. These hydrothermal cracks 
opened sufficiently due to fiberboard swelling to break the semi-ductile 
metallization on the cell front surface. Delamination cracks within the fiber 
board vyere also noted by edge examination of the hydrothermal ly aged cells. 

AC impedance measurements readily detect the damage effects of hydro- 
thermal aging and cell cracking as shown in the plots of Z' versus frequency 
in Fig, 23 and the Cole-Cole plots of Fig. 24. These curves show an approx- 
imate five fold reduction in the shunt resistance Rg}^ of the hydrothermal ly 
cracked solar cell. These abrupt and major changes in solar cell response are 
easily recognizable in a solar array as discussed in Appendix II, 

In order to understand the role that corrosion plays in a selected 
failure mechanism of Low Cost Solar Cells, surface analysis has been made 
using Auger electron spectroscopic profiling. A particular failure mode which 
was considered was that of cell cracking. Although cracking may or may not be 
induced by corrosive mechanisms, the resulting void produced by a crack allows 
a region for possible concentration of electrolytic impurities which aggravate 
the failure. Indeed, cracking of one cell may not necessarily produce an 
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Fig. 23 Log Z vs log w for EVA encapsulated cell before (|1) and after (#2) 
exposure to RH 100% at 75°C for 8 days. 
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Fig. 24 Cole-Cole plots of AC Impedance showing easily identified effects of 
silicon cell cracking. 





Rockwell International 

Science Center 


SC5106.104AR 


litimediate failure of an array, but the long term interaction of a corrosive 
environment could lead to a more serious defect at the crack by fonning of 
current carrying paths which would ultimately short circuit the cell. Work at 
JPL has shown that corrosion couples to the cracking process. 

In order to investigate this possibility, an Auger profile of a 
fracture surface adjacent to the metallization has been made. The sample was 
supplied by Dr. Gupta of The cell was cracked and had been exposed 

to UV light with intermittent soaks in distilled water. Figure 25 shows a 
schematic of the fractu>'e surface of the cracked cell. The front surface 
metallization of the solar cell consisted of a Pb/Sb solder whose major 
components were Pb, Sn and Fe with observable quantities of S, Cl, K, 0, Cr, 

Ni and C. The substrate (back surface) electrode is a Ni/Pb alloy. 

After a 200 sec sputter of superficial organic contamination, the 
fracture surface of the Si shows primarily Fe, Ni , and Cr distributed as 
illustrated in the inset in Fig. 25. The fracture surface of the Si also 
contained significant portions of C and traces of S, Cl, 0, and Cu. The 
distribution suggests a migration of Cr and Fe species from the upper elec- 
trode, possibly as a result of a corrosive mechanism. The sample also 
illustrates some migration of the Ni. The metals which have apparently 
migrated are expected to be the most active from a thennodynamic point of view 
since the thermodynamic tendency to electrochemical oxidation takes the 
following order Cr > Fe > Ni « Pb > Sn. although kinetic aspects of passi- 
vation and alloying influences upon chemical potential can produce significant 
changes in this ordering of actual lability, the selective migration of Fe and 
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Cr from the Pb/Sn electrode which exhibits positive polarity under illuiiin- 
ation argues for an electrochemical process leading to the dissolution of Fe 
nd Cr from the upper electrode. 

An additional HTSA computer run was conducted in which moisture 
exposure was suppressed at Oji R. H. and only the thermal cycle of lower 
Fig. 19 was imposed on the 4 layer solar cell model. The purely thermo- 
elastic stresses are shown by this computation to be noncritical and 
noncunul tive. 


5. 4 Interfacial Bonding for Corrosion Protection 

A good deal of evidence has been accunulated in the experimental 
studies of this program that environmentally resistant Interface bonding of 
encapsulant to all solar cell and metallization surfaces is essential for 
corrosion protection. This problen is even more acute when individual cells 
display microcracks as discussed in the proceeding sections. The atmospheric 
corrosion model discussed 1n earlier sections and summarized in Fig. 11 
defines explicit factors for reduction of corrosion rate by suppression of 
ionic diffusion mechanisms. 

In order to extend the theoretical range of this corrosion model, an 
additional acid-base criteria for interface bond stability, as shown in 
Fig. 26, will be explored. This criteria was originally developed and 
verified by Bolger and Krchaels^^^^ and has been incorporated into a general 
adhesion theory by Kaelble.^^^^ It is also evident that E. PI ueddemann^^^^ 
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ha'-( recently exploited this system of definition In selection and testing of 
silane adhesion promoters for LSA. The diagram of Fig. 26 provides a compact 
graphical method for selecting either addle (right margin) or basic (left 
margin) organic radicals to achieve respective high or with metal oxides 
(see upper margins) of differing isoelectric points for surfaces (TEPS). 
Interfaces formed according to the combined criteria Fig. 26 will be tested by 
atmospheric corrosion monitors and environmental corrosion simulation as 
already developed In the prior program. Additionally, a new ion migration 
test method, recently described by Mittal and Lussow^^®^ will be explored as a 
small scale screening test for Interface stability and corrosion resistance. 


53 

C3047A/ES 



SC80-9331 




Rockwell International 

Science Center 


O. 

CM 

X 


CO 

X 

o 

CM 

lo 

X X 


z 

o 

X 


O JP 

“O- 

• 0-^0 0 
o y o 
O § Q O 

o ? o o 

X O X X 


SC5106.104AR 


X 

M 

S 

o 

X 


II I III II 




no 


iv)Vj,d oiaiov 

o 


(Cl 


C*J 


CM 


(O 

I 


CM 

a.«» 

+* - 
a 
O 



CM 


03 


<0 


CM 


oisva 

Mil II I II 


X X N fM 
O Zi X 
^ II z 2: 
X " X ^ 

Z X X 

X o 

•©• 


2 «M 


CM WO 
^ X CM 
w Z X 

1 o 

2 o 

o X 

o 


CO 

z 

o 

h" 

O 

< 

X 

LU 

H 

Z 

u 

9 

o 

< 


< § 


<1 

>• 

H- 


X 
O 
X, 
a. ! 


< 

a 


CO 

a 

m 

II 

< 

<I 


X 

< 

H 

co 

a. 

O 

UJ 

(3 

Z 

< 

X 

w 


!fj 

V 

X 

a 

V 

< 

< 

a 


CO 

z 

o 

h“ 

U 

< 

X 

LU 

H- 


CO 

a 

CQ LU 

<3 

> 
h- 


X 

< 

X 

o 

X 

a 


X 

< 

a 


X 

< 


X 

< 

H 

CO 

a 

O 

X 

CP 

z 

< 

X 


X 

< 

3c: 

a 

V 

X 

a 

V 

X 

a 

LU 


54 


lEPS — 

ACIDIC BASIC 

$=pH— m 

^ 0.0591 

Fig. 26 Acid-base interface stability criteria. 





Rockwell International 


Science Center 

SC6106.104AR 


6.0 CONCLUSIONS 


Atmospheric corrosion monitors (ACM) have been returned to the 
Science Center following 13 months deployment at the Mead, Nebraska test 
site. During this period the ACM units recorded the corrosion protection 
function of an encapsulant system and correlated Mead cl imatol ogy with 
corrosion rates of a noriencapsulated ACM. The fundamental assumptions of a 
new atmospheric corrosion model were verified in this study. This corrosion 
model predicts that corrosion rate is the product of a condensation proba- 
bility (P^) and the maximum ionic diffusion current Encapsulant 

corrosion protection is specifically related to its efficiency in suppresssing 
II at the potential corrosion interface. 

AC impedance measurements combined with impedance spectriiii analysis 
and feed back control to a branch circuit appear as a direct means of remotely 
locating degraded LSA modules and modifying series to parallel S/P inter- 
connects between modules to achieve maximum LSA power efficiency. 

A hydrothennal stress analysis (HTSA) computer model has been 
successfully applied to evaluate the combined effects of temperature- humidity 
cycling on development of internal stresses and solar cell cracking in solar 
cell modules using a fiberboard substrate. The effect of protective polymer 
coatings of ethylene vinylacetate (EVA) and. polytrifl uoro-chloroethylene 
(Kel-F) is to delay but not change the failure process. 
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A new materials selection criterion for encapsulant bonding is being 
developed which incorporates both ionic cor)duction and el ectrochetnical 
mechanisms of corrosion. Tlie new selection criteria and test methodology are 
specifically directed at suppressing the micro-corrosion process in the 
presence of internal defects such as micro-cracks. 
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7.0 RECOMMENDATIONS 


The following specific recommendations In the area of corrosion and 
encapsul ant life prediction are presented from this study: 

1. AC Impedance monitoring combined with spectrum analysis and feed 
back control is suggested for field trails at an ISA test site. 

2. Develop materials selection criteria and tests for corrosion and 
environment resistant interfaces with 20-year life capability. 

3. Develop and validate a more general reliability and life 
prediction model for ISA which Includes corrosion, Interfacial 
integrity and hydrothermal stress analysis (HTSA) as specialized 
subjects and directly describes photovoltaic current- vol tage 
(I-V) response. 
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8.0 NEW TECHNOLOGY 

This study has developed arid demonstrated new corrosion models and 
test methods. AC impedance measurements and analysis has been advanced to a 
potential field deployable state for nondestructive evaluation (NDE) and 
remote manageiient of ISA branch circuits. Hydrothermal stress analysis (HTSA) 
has been demonstrated as a method for preilicting solar cell life time and 
failure mechanisms. Atmospheric corrosion monitors (ACM) are effective tools 
for evaluating atmospheric corrosion resistance of encapsulant bonding 
materials and evaluating climatology influence on corrosion rates. 
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APPENDIX I 

MATHEMATICAL RELATIONS FOR CURRENT -VOLTAGE (I-V) RESPONSE 

IN SINGLE SOLAR CELLS 


Introduction 

It is often useful to develop an analog model of physical response in 
order that the interrelations between variables are more fully understood. 

Tlie model for I-V response presented here is mathematically simple and yet 
resembles many aspects of real physical response of single solar cells. Of 
particular interest in this discussion is the analytical description of the 
fill factor (F), and the maximun power point of the I-V curve which 
describes cell power generating efficiency. Table 1 summarizes the symbols 
and nomenclature used in this discussion. 

General Relations 

Let us assume that short circuit current Iq is proportional to 
insolation S: 


= Id + KjS 


( 1 ) 


and that the open circuit voltage Vq decreases with increased temperature T 
and insolation S as follows: 
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l<2 


( 2 ) 


where Kj, Kg, and K 3 are positive constants. By introduction of a shape 
factor K, which is defined as follows: 


K 



An 


V^. 


(3) 


where R 3 and are respective series and shunt resistance the IV curve can 
be expressed as follows: 


I = I(, [1 - exp K(V - V^j)] 

(4) 

X = 1 - exp{ - KVjj) 

(5) 


Determining Rg and R ^^ 

By differentiation of Eq. (4) for constant Iq, V^, and K we obtain 
the following relation: 

At the open circuit condition, where V = Vq and I = 0, Eq. ( 6 ) becomes: 

(7) 
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At the short circuit condition, where V » 0 and I * Io> (6) becomes 

^ i:exp(- KVg)] » (8) 

The limiting slopes of the I-V curve are defined by Eq. (7) to determine Rs at 
V = Vq and Rg + R 5 ^ at I * Iq. 


fletennining the IV Fill Factor F 

The fill factor of an IV curve is defined by the following integral 

rel ation; 


F = -jV / I ‘IV 
0 0 0 


(9) 


By substituting Eq. (4) into Eq. (9) we obtain; 


F =^/ [1 - exp K(V - V )] dV 
0 0 

-1 ^0 ^0 

<1V -// Cexp K(V - V )] dV 
0 0 0 0 


( 10 ) 


Integrating Eq. (10) by parts we obtain the following relation: 


F = Cl - exp(- KV^)] 


( 11 ) 


Substituting Eq. (5) into Eq. (11) provides the following relation: 
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p , ^ ^2) 

0 

as shown by Eq. (12) when KV^ « 0 the fill factor F « 0,50 and when KVq 
achieves high values KV^ > 50 the fill factor approaches unity. 


Determining the Maximum Power Voltage 

The power output W of the solar cell is determined by the Following 

product: 


W = IV (13) 

The maximum power point, with respect to voltage V, is obtained by 
differentiating Eq. (13) as follows: 

d W _ T 1 If d I / , . \ 

w ~ ^ ^ W 

By solving for maximum power where dW/dV = 0 and in substitution of Eq. (4) 
and Eq. (6) into Eq. (14) we obtain the following expression: 

HU VK 

^ = ^ [1 - exp K(V - Vq)] - K exp K (V - V^) = 0 (15) 

By simplification and rearrangement, Eq. (15) defines the voltage where 
maximum power VI, ^ is produced from the I-V curve by the following relation; 
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Xn (KV^ + 1) KVj„ » KVq (16) 

Solutions to Eq. (16) are obtained by introducing values of KV^, in the left 
side and solving for KVq. The subsequent determination of ratio of V^j/Vo 
versus KVq is obtained by dividing the Input on the left side of Eq. (16) 
by the output KVq on the right side of Eq. (16) to obtain the dimensionless 
ratio VVq. 


Determining the Maximum Power Output 

For a given I-V curve the maximum power output is expressed by the 

rel ation: 


W “ I V_ 
m mm 


(17) 


where and are the respective current and voltage at the maximum power 
point. The relation for 1,^ is provided by substitution of V = Vp, in Eq. (4) 
to provide the following relation; 

[1 - exp K(V|„ - Vj,)] (18) 

By substituting Eq. (18) Into Eq. (17) provides the following relation for 
inaximuTi power: 

\ - “P <19) 
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Alternatively, the maxlinun power can be reduced by the product IqVq to 
define a maxinuin power output ratio as follows 

W V 

^ “ TV- • -? Cl ■ “P - VJ] (20) 

"0 0 0 

Shunt to Series Resistance Ratio 

By rearranging Eq» (3) the ratio of shunt to series resistance is 
given by the following relation: 

^•exp(KV^)-l (21) 

s 

Sample Computations for IV Properties 

The procedure for conducting abbreviated calculations for IV response 
proceeds as shown in Table 2. The calculation proceeds from left to right 
across each line of Table 2 and from top to bottom through the list of 
equations in lower Table 2. For convenience, this calculation simply inputs 
an assigned value of KVp, into Eq. (16) and computes a value of KVq in the tw 
left columns of Table 2. The dimensionless ratio Vj^^/Vq is defined and the 
reciprocal value calculated as shown in the third and fourth columns of 
Table 2. The remaining properties F, P^, and R$h/Rs in the right column of 
Table 2 are readily calculated to complete the calculations. 

It is evident that in full calculations both Eq. (1) and Eq. (2) 
would be utilized with input of two environment conditions T and S and six 
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solar cell constants Ip, Vp, K» K^, Kg, K 3 with computation of Ig and Vg. The 
next computation inputs given value of KVg into Eq. (16) and successive values 
of substituted to satisfy the equality with a determined value of as 
output. Successive steps of the full calculation follow those of the 
abbreviated procedure of Table 2. 


Table 1 

List of Principal Current-Voltage (I-V) Design and 
Performance Parameters 


Symbol 

Meaning 


Design Parameters 

Coefficient of light developed current 

^2 

Thennal coefficient of voltage 

K3 

Thermal coefficient of light 

'0 

Dark current 

Vd 

Dark vol tage 

Rs 

Series resistance 

^SH 

Shunt resistance 

T 

Environment Parameters 
Ambient temperature 

S 

Ambient insolation 

^0 

Perfonnance Variables 
Short circuit current 

Vo 

Open circuit voltage 

L, X 

Current- vol tage (I-V) response factors 

F 

I-V fill factor 


Voltage at the maximum power point 


Wattage at the maximum power point 

^m 

Current at the maximun power point 

Pm 

Proportion of ideal power output 
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Table 2 

Sample Computations of Solar Cell IV Response 


KVm 

KVo 

Vra/Vo 


F 

'^m 


.001 

1.9995E-3 

.5001 

5.006E2 

.50017 

.2501 

.002 

.002 

3. 9978E-3 

.5003 

2. 506E2 

.4991 

.2503 

.004 

.005 

9. 9875E-3 

. 5006 

1.006E2 

.5007 

.2506 

.0100 

.01 

1.995E-2 

.5012 

5.0627E1 

.5016 

.2512 

.0201 

.02 

3. 903E-2 

.5025 

2.6125E1 

.5033 

.2525 

.0406 

.05 

9. 789E-2 

.5061 

1. 0724E1 

.5082 

.2562 

.1038 

.10 

1.953E-1 

.5120 

5.6366 

.5162 

.2623 

.2157 

.20 

3.823E-1 

.5231 

3. 1474 

.5318 

.2744 

.4657 

.50 

9.055E-1 

.5522 

1.6789 

.5744 

.3090 

1.473 

1.00 

1.693 

.5906 

1. 2254 

.6347 

.3618 

4.436 

2.00 

3. 099 

,6455 

1,0472 

.7245 

.4506 

2.117E1 

5.00 

6.792 

.7362 

1.0011 

.8539 

.6142 

8. 895E2 

10 

1.2398E1 

.8066 

1.0000 

.9193 

.7333 

2.423E5 

20 

2.3045E1 

.8679 

1. 0000 

.9566 

.8266 

1.019E10 

50 

5.3932E1 

.9271 

1.0000 

.9815 

.9089 

2.6442E23 

100 

1.0462E2 

.9559 

1. 0000 

.9904 

. 9465 

2. 728E45 

200 

2.0530E2 

.9742 

1.0000 

.9951 

.9693 

00 

500 

5.0622E2 

.9877 

1. 0000 

.9980 

.9857 

00 

1000 

1.00691E3 

.9932 

1.0000 

.999 

.9932 

00 


KVq = iln (KVj„+ 1) + 

= 1/[1 - exp(-KVQ)] 

F = - (KV )-l 

Pm= (Vri/Vq)^"^ Cl - exp K(V^- V^l)] 

R^h/Rs "" ^ 


IV Shape Factor (Eq. 16) 

Max. Power Voltage Ratio 
Function of Shape Factor (Eq. 5) 
IV Fill Factor (Eq. 12) 

Max. Power Ratio (Eq. 20) 

Shunt/ Series Resistance (Eq. 21) 
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APPENDIX II 

Analysis for Frequency Dependence of AC Impedance 


Introduction 

AC (alternating current) Impedance measurements have been applied exper- 
imentally to ISA (low cost solar array) single solar cells or to modules of 50 to 
75 cells electrically connected in series. The AC impedance measurement is 
conducted by scanning the frequency spectriin from 1 Hz through 1 MHz using low 
amplitude AC voltage (±10 nV) and constant temperature. The frequency response 
for AC impedance of a single solar cell appears to fit the classical Debye type^ 
single relaxation mechanism. When multiple solar cells are examined, the 
frequency response for AC impedance is more complex but appears to be well 
described by a symmetrically broadened frequency distribution of Debye type 
relaxation processes. The empirical distribution function developed by Cole and 
Cole provides a set of general equations for broadened AC impedance response 
which reduce to the Debye single relaxation as a special case. In this 
discussion, the analysis centers on the Cole-Cole distribution function as a 
means for evaluating performance parameters of both single solar cells and 
multiple cells in modules or branch circuits. The AC impedance measurement 
represents a readily conducted nondestructive test which when combined with the 
dielectric analysis discussed here becomes a method for conducting nondestructive 
evaluation (NDE) of ISA performance. A list of the principal parameters of this 
analysis is presented in Table I. 
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Tlie Cole-Cole Distribution Function 

The characteristic relaxation time Tq for three element circuit 
‘^4io\vn in Fig. 5 is defined by the following relation: 


Tq » 



( 1 ) 


The storage or in-phase impedance z' is defined as a function of test frequency u) 
(radians/s) by the following relation: 




Rsh cos(3ir/2)] 


1 + (oit) ^ + E(wt)^ (cos 


(3 


( 2 ) 


The loss or out-of- phase impedance z" is defined as a function of test frequency 
by the following relation: 


z 


II 


1 + {0)xr^ 


sin(|3ir/2) 

+ 2(wt)^ cos($ir/2) 


(3) 


The distribution quality factor 0 is unity, 6 = 1.0, for the Debye single 
relaxation and defines a range 0<$«sl.0 for symmetrically broadened distribution 
with multiple relaxation times. 

Tlie Cole-Cole distribution function <i>(ln t) , is defined by the following 

rel at1on: 
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^(Inr). (4) 

coshCe In (—)] + cos Bn 

The distribution function Eq. (4) fulfills the following general normalizing 
rel ation: 


I S(ln t) 


d In T = 1 


The detailed form of the distribution function of Eq. (4) relates to the dis- 
tribution quality factor B as will be explored in subsequent discussion. 

Combining Eq. (2) and Eq. (3) to eliminate ut provides the following 
equation for a circle: 

2 

( 6 ) 


..("s * 

2 

X 

R-u cotan (Bn/2) 

2 

Rgl^ cosec Bn/2 


T 

^ ^ 2 


2 


The radius r^. of this circle is: 


/Bttv 

r^ - 2 cosec (2 ) 


(7) 


and its ce ter has coordinates: 


r' = 






. Y cotan (f ) 


( 8 ) 

(9) 
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AC impedance data that conform to the Cole-Cole distribution function can be 
conventionally plotted as a circular function on z' and z" coordinates. 

Distribution Analysis of AC Impedance Data 

A nuiierical analysis of AC impedance has been developed by Messrs. Gene 
Meyer and Gary Lindberg of the Science Center for a least squares circle fit of 
the Cole^Cole distribution. A schematic representation of the circular least 
squares fit is illustrated in Fig. 9 where r^,, z^, and z'^ are respectively 
defined by Eqs. (7-9). 

The computation is addressed to reducing the distance d^ between the 
best fit circle aixl the experimental points represented by zl , z!j. From Fig. 9 
v/e have: 


“ tXi' 




( 10 ) 


*1 ° K ' ’’c 

yj = Z'j - zj; - sin 0. (12) 

where Xj and y^* are the respective z' and z" axis errors between the curve and 
the point z*. , z'j. To obtain the best fit circle to a set of AC impedance data 
points the objective is to minimize D which is defined from Eq. (1-3) as follows; 
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" 2 
D " E d/ 

Ul ^ 


’ • ’'i • •'c 


(zlj ~ - »'c 


By taking the partial deviatives of D in Eq. (13) with respect to xj,, x"^ and r^, 
and setting 3D/3x' = 0, 3D/3x", and 3D/3r. equal to zero the following relations 

C w w 

are obtained; 






n 

E cos 0. 
i=l ^ 


n 

E z!. = 0 
i=l ^ 


(14) 


q = 



n 

nx^ + r^ E sin 0. 
c c 1 


n 

E zli' = 0 
i=l ^ 


(15) 


h = 



n n 

nr^ + z' E cos 0. + z" e sin 0. - 
c c 1 c 1 


n n 

E z'. cos 0. - E z!;' sin 0 . = 0 (16) 

i=l ^ ^ i=l ^ ^ 


TJ L"--. 
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Inspection of Fig. 9 periods the following relations for sin in Eqs. (13-16), 


sin = 




c(2; - z;)'' + (z'j - 


(in 


cos 0^. 


z'l - z" 
1 c 


Cz; - z^,)‘+ (z'j 


(18) 


Inspection of Eqs. (14-16) shows that: 





n 

t cos 0 . 
i=l ^ 



i=l 


sin 0. 



^ 0 h 

2 cos 0*; 

i=i ' 


n 

Z sin 0z; 
i=l ^ 



n 


A set of three equations and the unknowns is expressed by Eqs. (14-16). Solu- 
tion of these equations can be acconipl i shed by use of a three parameter analog 
of the Newton-Raphson method as described by Henrici.^ Expanding the example 
of two equations and two unknowns described by Henrici to our case, provides 
the following determinant relations, where K and K + 1 represent successive 
iterations where xj. = ^c ^^c^K ^'c 
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3f 

3f 

3f 

axT 

3F 

3r_ 

c 

c 

c 


air 

i9_ 

3r_ 

c 

c 

c 

3h 

3h 

3h 





(19) 



+ 1 “ 








3h_ 


( 20 ) 


8f 


c^K + 1 




3h 


-f 


-h 


8f 

ar. 


ia_ _n ia_ 

^ y 9r 

c c 


3h 

ar. 


M. 


( 21 ) 
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(rJ 


c^K + 1 




8f 

3f 

3x1 

IxT 

c 

c 


% 

c 

c 

3h 

3h 

aTr 

c 



-f 


-h 


( 22 ) 


The computation initiates at K * 1 by selection of tl'u'ee widely separated 
experimental points (z'., z^ in the graph of Fig. 9 in order to make a good 
initial estimate of xj., xjl and r^ for K - 1 by triangulation. 


Using these values of xj., x^, and r^. one sums ove* the data set to 
obtain the following s unmat ions: 


which remain constant through subsequent iterations. Additionally suiimations 
are made for each iteration K for the following: 


n n n n 

E sin 6.; Z cos 6.; Z z'. cos 9,.; Z z'l sin 0. 

i=l ^ i=l ^ i=l ^ ^ i=l ^ ^ 

which will correlate with the introduced values of (z')i,, (z")t,, {r)v into 

C l\ C K C In 

Eqs. (17-18). The above sunmations are introduced into Eq. (14-16) to solve 
for the K iteration values of f, g, and h and M|^ from Eq. (19). 
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The above K iteration coefficients are applied in Eqs. (19-22) to 

obtain the K + 1 iteration solutions for x' , x" and r_. These iterative cal- 

c c c 

culations continue until terminated by an appropriate convergence test condi- 
tion, such as;^ 


0 < 


^ " ‘^M - l 


< e 


K+1 


(23) 


Where Q represents the respective coefficients z' or z" or r. for the K+1, K-, 

c c c 

and K-1 iterations and the assigned convergence condition for obtaining a 
final solution where: 


Q = (24) 

and Q is the best value of the respective coefficients z^, z^ and r^. The 
iteration calculation is completed by satisfying Eq. (23). The radial stan- 
dard deviation 5^, of experimental points to the calculated radius r^ is 
obtained by utilizing the parameter D from Eq. (13) in the following conven- 
tional relation 


Summary 

The above definitions and computational procedures provide a means of 
isolating significant solar cell performance parameters. For example, by mea- 
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suring a limited frequency range of AC Impedance response this analysis per- 
mits direct calculation of the shunt and series Rj resistances of the ISA 
module or branch circuit. By direct measurement Rg and R£^ are evaluated as 
the respective high and low frequency asymptotes of a broad range of frequency 
measurements. This analysis also provides a direct method for evaluating the 
unifonnity of single cell responses within the solar array in terms of the 
Cole-Cole distribution function. This analysis is currently incorporated in a 
micro- computer program for AC impedance evaluation of LSA. 
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Table 1 

List of Principal AC Impedance Measurements 
and Perfonnance Parameters 

Symbol Meaning 

Measurement Parameters 

2 ' Storage component of AC impedance at frequency w 

z" Loss component of AC impedance at frequency w 

w = Zrtf Test frequency (rad/s) 



Impedance Respon se 

Shunt resistance 
Series resistance 
Capacitance 

Cole-Cole shape factor (0<3«;1.0) 
Characteristic relaxation time 

Parameters of Cole-Cole Numeric Analysis 

Mean radius of Cole-Cole arc 
Center of storage axis of Cole-Cole arc 
Center of loss axis of Cole-Cole arc 
Standard deviation of the mean radius r^ 
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